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Abstract 

We make a perturbative calculation of neutrino scattering and absorption in hot and dense hy- 
peronic neutron-star matter in the presence of a strong magnetic field. We find that the absorption 
cross-sections show a remarkable angular dependence in that the neutrino absorption strength is 
reduced in a direction parallel to the magnetic field and enhanced in the opposite direction. This 
asymmetry in the neutrino absorbtion can be as much as 2.2 % of the entire neutrino momentum 
for an interior magnetic field of ~ 2 x 10 17 G. We estimate the pulsar kick velocities associated 
with this asymmetry in a fully relativistic mean-field theory formulation. We show that the kick 
velocities calculated here are comparable to observed pulsar velocities. 

PACS numbers: 25.30.Pt,21.65.Cd,24.10.Jv,95.85.Sz,97.60.Jd, 
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I. INTRODUCTION 



Hot and dense hadronic matter is a topic of considerable current interest in nuclear 
and particle physics as well as astrophysics because of its associated exotic phenomena. In 
particular, many studies have addressed the possible exotic phases of high density matter. 
Neutron stars are thought to be among the most realistic possible sites for the presence of 
such high density matter. For example,the possible existence of an anti-kaon condensation 
in neutron stars has been suggested [lL and the possible implications for astrophysical 



phenomena have been widely discussed 2 



In this context, a number of works have discussed effects of a softeningo 



: the equation of 



3§. Other work 

14 3- 

Among 



state (EOS) on both the static and dynamic properties of neutron stars la- 
has considered the thermal evolution of neutron stars by neutrino emission 
them, Reddy et al. [20[ have studied neutrino propagation in proto-neutron stars (PNSs) as 
a means to examine the hyperon phase in the high density region. 

On the other hand, magnetic fields are thought to play an important role in many astro- 
physical phenomena such as as ym metry in supernova (SN) remnants, pulsar kicks 2l|, and 
the existence of magnetars [22j, [23|. Indeed, strong magnetic fields seem to be a crucial part 
of the still unknown mechanism for non-spherical SN explosions and the origin of the high 
velocity [24J that proto-neutron stars (PNS) receive at birth. 

Although several post-collapse instabilities have been studied as a possible source of 
non-spherical explosions and pulsar kicks, the unknown origin of the initial asymmetric per- 
turbations and the uncertainties in the numerical simulations make this possibility difficult 



to unambiguously verify 



25 



26| . Another viable candidate is the possibility of asymmetric 



27 



28j or as a 



neutrino emission either as a result of parity violation in the weak interaction 
result of an asymmetric magnetic field [29] in strongly magnetized PNSs. Recent theoretical 
calculations 30, |31[ have suggested that even a ~1% asymmetry in neutrino emission out of 
the total neutrino luminosity of ~ 10 53 ergs could be enough to explain the observed pulsar 



kick velocities. 

Reddy et al. [2^ studied neutrino propagation in PNS matter including hyperons. How- 
ever, these models did not include the effects of a strong magnetic field. Although Lai et 



al. 



30. 



3l| studied the effects of a magnetic field on asymmetric neutrino emission, the 



neutrino-nucleon scattering processes were calculated in a non-relativistic framework 30]. 
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In Ref. 



32] we reported for the first time our calculated results of the absorption cross- 



sections in hot dense magnetized neutron-star matter including hyperons in fully relativistic 
mean field (RMF) theory ^i,]M\. 

In this paper, we give more detailed information of the neutrino scattering and absorption 
cross-sections in magnetized neutron-star matter studied in the RMF theory. In the present 
RMF framework we take into account the Fermi motion of the baryons and electrons, their 
recoil effects, the distortion effects of Fermi spheres made by the magnetic field, and the 
effects of the energy difference of the mean field between the initial and final baryons. 
We then solve the Boltzmann equation for neutrino transport in a ID model and discuss 
implications of the present results for the pulsar kicks. 

In Sec. II we introduce the EOS of nuclear matter based on the RMF theory. In Sec. 
Ill we explain the neutrino scattering and absorption cross-sections in baryonic matter with 
the strong magnetic fields. Then, we show numerical results for neutrino reactions and 
propagation in matter at finite temperature and discuss the possible associated pulsar kicks 
of magnetized PNSs in Sec. IV. In Sec. V we summarize our work in this paper. 
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II. NEUTRON-STAR MATTER IN THE RELATIVISTIC MEAN-FIELD AP- 
PROACH 



First consider neutron-star matter including nucleons, Lambdas, electrons and electro- 
neutrinos (v e ) without a magnetic field in the RMF approach. The Lagrangian density for 
the matter field is written as 

£ Matter = &RMF + £-Lep, (1) 

where the first and second terms are the RMF and lepton parts, respectively. The lepton 
part includes only the kinetic energy of the electrons and the electron neutrino v e . The RMF 
part is given by 

Crmf = i'Niil^ - M N )ip N + g^NipNO' + Quj^nI^n^ 
~ 1 C 

- U W\ + 2 m l UJ ^' 1 ~ -^2^Nl^T a llj N ^ N ^T a llj N} (2) 

where ipN, i>A, an d u> are the nucleon, Lambda, sigma-meson and omega-meson fields, 
respectively, M^r, M\, and m u are the masses of corresponding particles, 7 M is the Dirac 
matrix, and r a is the Pauli matrix in the isospin spaces. U[a] is the self-energy potential of 
the scalar mean-field given in Refs. Q,Q. The last term describes the vector-type iso- vector 



interaction between two nucleons and is equivalent to p-meson exchange 33j. [Note that 
here and in the following we adopt natural units, i.e. H — c—1] 

From the Euler-Lagrange equations of the above Lagrangian, the Dirac spinor of the 
baryon Uf,(p, s) is obtained as a solution to the following equation: 

li>-MZ-Uo(b) lo }u b (p,s)=0, (3) 

where Ml is the baryon effective mass, and Uo(b) is the time component of the mean-field 
vector potential. We hereafter introduce the Feynman dagger p 7 = 7 M p M for convenience. In 
the RMF approach the nucleon and Lambda effective masses are given by 

M* N = M N -U„(N), 

Ml = M A -U S (A), (4) 
with the scalar mean-field potentials: 

U s (N)=g a (a), U.(A) (5) 



The scalar mean-field (a) is given by 

d 



U[(a)} = Qcr [p s (p) + Ps(n)] + g^p s (A), (6) 
with the scalar densities, p s (b) {b = p,n,A), defined by 

Ps(b) ^^yj^P [^V^ip)] + nt\^\p)]} (7) 

where e£ are the single particle (+) and antiparticle (— ) energies, E£(p) = \fp 2 + M fe * 2 , 
and the Fermi distributions, ra[ ± ' ) (e[ ± ' ) ), are defined as 

= L ( § ) 

l + exp[( e W± £6 )/T] 

with the temperature T and the chemical potential e^. In addition, the baryon single-particle 
energy is written as e ( f z \p) = E£(p) ± U (b), and the U (b) are given by 

U (p) = + Pn)+gtPh} + inPi{Pp- Pn), (9) 

m ui 1V1 N 

U (n) = % {g u (p P + Pn) + g£,Pk} - tS"(Pp- Pn), (10) 

m uj IV1 N 
U (A) = { 9u} ( Pp + p n ) + g£p A } (11) 

in terms of the proton, neutron and Lambda number densities, p p , p n and p\. 

In this work, neutron-star matter at finite temperature will include protons, neutrons, 
Lambdas, electrons and neutrinos with the conditions of charge neutrality and beta equi- 
librium being satisfied. Then, the proton number density is equal to the electron number 
density, p p = p e , and the chemical potentials obey the following condition 

£n — £ a — £ p + e e • (12) 

In addition, we assume that the lepton fraction is fixed to be Y L = (p e + p u )/ P B, where p v 
is the neutrino number density, and ps = p p + p n + Pa- 

In the actual calculations we use the parameter sets of PM1-L1 [35(], which give the binding 
energy per baryon BE = 16MeV, the effective mass M^/M^ = 0.7 and the incompressibility 
K = 200MeV at po= 0.17 fm~ 3 in nuclear matter, and the sigma- and omega-Lambda 
couplings are 2/3 of those for the nucleon, = 2/3q u,,,- A similar relation is found 



naturally within the quark meson coupling (QMC) model 40 ]. 
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FIG. 1: (Color online) Upper panels (a) and (c) show the density dependence of the total energy per 
baryon Et/A in neutron-star matter for T = 20 MeV (a) and 40 MeV (c). Solid and long-dashed 
lines represent the calculated results in systems with and without Lambda particles. The lower 
panels (b) and (d) show the number fractions of protons x p , Lambda particles x\, and neutrinos 
x v for T = 20 MeV (b) and 40 MeV (d). Solid, dot-long-dashed, and short-dashed lines represent 
the calculated proton, Lambda, and neutrino number fractions, respectively, in a system including 
p, n and A. Long-dashed and dotted lines represent the calculated proton and neutrino number 
fractions in a system without Lambdas. In the present calculations we use the parameter-set 
PM1-L1 [35| for the RMF and the lepton fraction is set to Y L = 0.4. 



In Fig. [T]we show the energy per nucleon in the upper panels (a and c) and the proton and 
Lambda fractions in the lower panels (b and d) at T = 20MeV (a and b) and T = 40MeV 
(c and d). In these calculations the lepton fraction is taken to be Yl = 0.4. The solid and 
dashed lines represent the results for matter with and without Lambdas, respectively. In 
addition the dot-dashed lines in the lower panels indicate the Lambda fraction. We see that 
the Lambda fraction appears when ps > 2p and affects remarkably the EOS for p# > 3po- 
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III. CROSS-SECTIONS FOR NEUTRINO REACTIONS IN MAGNETIZED 
PROTO-NEUTRON STAR MATTER 



We separate the Lagrangian density into the following parts 

£ = &RMF + £-Lep + £-Mag + £-W ■, (13) 

where the first, second, third and fourth terms are the RMF, lepton, magnetic, and weak 
interaction parts, respectively. Detailed expressions for these terms are explained in the 
following subsections. 



A. Dirac Wave Function in a Magnetic Field 

We assume that there is a uniform magnetic field along the z-direction B = Bz and 
for £?<10 18 G, the effect of the magnetic field on baryons is small enough to be treated 
perturbatively. The magnetic part of the Lagrangian density is written as 

£-Ma g = &BM + ^eMi (14) 

where the first and second terms describe the magnetic interactions of baryons and electrons, 
respectively. 

Ignoring the contribution from the convection current and considering only the spin- 
interaction term, the baryon magnetic-interaction Lagrangian density can be written as 

C-bm = ^UbiibCr^ipbF^ = ^2/j, b i> b a z ^ b B (15) 
b b 

with the electromagnetic tensor given by F^ u = d^A u — d u A 1 *, where is the electro- 
magnetic vector potential, a^ u = [7^,7^/2?, <j z = diag(l, — 1, 1, — 1) and //& is the baryon 
magnetic moment. We then obtain the baryon wave functions by solving the following Dirac 
equation: 

$ - M* - U (b) l0 - fi b B lo a z ] u b (p, s) = . (16) 

The single particle energies e b (p, s) and the Dirac spinors in the limit of a weak magnetic 
field are given as follows: 

e b (p,s) = sjpl+t^JpY+Mf + ^B^j +U (b) 

» E* b (p) + U (b) + AE;(p)s , (17) 



with 

\ VP) 



and 



u b (p,s)u b (p,s) = —^-[E* b (p) lo -p- 7 + M*](l + s^(p)) 



+ ^ vrr=m ("^ (p)7 ° + M + ^ 3 - p ^ ~ p ^ ' (19) 

8^*(p)VPT + M 

with 

a(p) = (a , or, a 2 ) = — ===^(p„ 0, 0, £ 6 *(p)) . (20) 

Detailed derivations of these expressions are presented in the appendix [A] The second and 
third terms of Eq. ( fl9l) do not appear in the non-relativistic framework, but their contribu- 
tions are negligibly small and omitted in the present work. 

When we calculate the electron contribution in Eq. (TH|) . we have to use another treat- 
ment. This is because electron mass is very small, and its current is almost a Dirac current: 

C eM = -<$ e 7 M ^M", (21) 

where ip e is the electron field. Also, the effect of a strong magnetic field on electrons may 
not be a small perturbation. The electron energy in the presence of a strong magnetic field 
is generally given by 

e e {n,k z ;s) = ^/k 2 z + m 2 e + eB(2n + 1 + s), (22) 

where n enumerates the Landau levels of the electrons. 

In the limit of B — > 0, the electron wave function becomes a plane wave. Therefore, we 
can use the same expression as Eq. ( f!9|) for electrons, except for the spin vector. The upper 
component of the electron Dirac spinor is an eigenvector of the matrix a z . The spin vector 
in the rest frame of the electron is then (0; 0, 0, 1). In the matter frame the boosted spin 
vector can be written as 

a(k) = a e (k) = (— , — r, 1 + — tttt^ , (23) 

K ' K ' W m e (E e (k)+m e y m e (E e (k) + m e ) ) 1 ; 

where k z and fcr are the components along the z-direction and perpendicular to the z- 

direction, respectively. 
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When \JleB <C e e , the summation over n can be approximated as an integration over 
energy, i.e. 

^ ~* ^l dXT ' ^T = 2eB{n+ 1 -)) . (24) 

n 

Note that the variable xt corresponds to p\ in the limit of B — >■ 0. Then, the expectation 
value of a operator O is given by 

2eB 



<0>- 



f9 i2 ^E / dk zn e (e e (n,k z ,s))(D(n,k z ,s) 

^ ' s n 

~ (2^)2 y rf;rT y d kzn e (e e (xT,k z ,s))0(x T ,k z ,s) 

where the electron energy is approximately e e ~ \Jk 2 + m? e — eBs. 

Actual calculations are performed in the limit of m e — > 0, so that the electron energy and 
the spin vector can be approximated by 

/To n e Bs . , , m e „ . . 

e e ^A/fc 2 + m2 — « fe + 7T>e£s , (26) 

v 2A/fc 2 + m2 |«| 

^^('-whJ' (27) 

where /x e = —e/2m e . 

In this paper the temperature is taken to be rather low relative to the baryon rest mass, 
T < 40 MeV. In this case the anti-lepton and anti-baryon contributions are negligibly small 
in the neutrino reactions. Therefore, we ignore the contributions from antiparticles, and 
omit the superscript '+' in the single particle energies e^\p) and the Fermi distribution 
n ( t\p, s). 



B. Neutrino Reaction Cross-Sections 



We consider neutrino reactions in neutron-star matter consisting of electrons and baryons 
(i.e. protons, neutrons and Lambdas). We assume individual collisions between the initial 
neutrino and the constituent particles, and calculate the neutrino scattering [y e — > u e ) 
and absorption {y e — > e~) cross-sections in the impulse approximation. Furthermore we 
only consider rather low temperatures, T <C as discussed above, and then ignore the 

9 



contribution from antiparticles. In addition, we treat this system as partially spin-polarized 
owing to the magnetic field. The cross-section can then be described in terms of the initial 
and final lepton momenta as fej and kf, as 

d 3 a Gp v-^ v-^ f d 3 pi 



n 4~ Ci C • Br J \ J 



dk 3 f 16n 2 

x rnieiipi, Si))[l - n f (e f (p f , s f ))} 

x (2n)5 3 (k i + Pi - k f - p f )S{\ki\ + Ci - e t - e f ), (28) 

where V is the volume of the system, and the index / denotes the final lepton species. The 
indices i and / denote the initial and final particles for both baryons and electrons which 
comprise the neutron-star matter. The function Wbl in Eq. ( 1281) is defined as 

W BL = J . LTNp, (29) 

4N|fc/|£*(P;)£/(P/) 

with 

LT = ^Tr{(^ + mO(l + 75M)r(l -75)^(1 -75)} , (30) 

and 

iV = -Tr{(j5 / + M / *)(1 + js4fSf)%(c V -c A j 5 ) 

x (jj. + M*)(l + 75^Si)7,(cy - c Alb )} , (31) 

where mi is the mass of the final lepton. 

Since we assume that the magnetic field is weak, the Fermi distribution and the delta 
function in the above equations can be expanded in terms of the magnetic field B. Then the 
cross-section can be presented as a sum of two contributions with and without the magnetic 
field 



d 3 a s , A _ <P°%A , d 3 Aa s , 



dk j dk j dk ^ 

where the indices S and A indicate the cross- sections for scattering and absorption, respec- 
tively. Note that a° s A is independent of B, and A.as,A is proportional to B. By considering 
the spin-dependence, we can express the Wbl as follows 

W BL = W + W iSi + W f s f + W e s t + W 2 s iSf + W 3 s lSi + W A s lSf . (33) 

Note that W e , W3 and W4 only appear when the final lepton is an electron. 
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When \nbB\ <C e b — Uo(b), the baryon Fermi distribution function can be expanded as 



n b (e b (p, s)) « n b (E* b (p) + U (b)) + n' b (E* b (p) + U (b))AE b (p)s, (34) 

and the electron distribution is written as 

n e (e e (k)) « n e (|fe|) + n^fe)^^!, (35) 

|fc| 

where ^(x) = dn b (x)/dx. In addition, the energy delta-function is also expanded as 

5(\ki\ +e i {p i ,s i ) -ei(k f ,si) -e f (p f ,s f )) 
« 5(N + + l/ (i) - 1%! - E}(p f ) - U (f)) 

+6'{\ki\ + E*{ Pi ) + C/ (z) - \k f \ - E}(p f ) - U (f))AE , (36) 

where 5 f (x) = dS(x)/dx, and 

771 

= AE^p^Si — AEf(pf)sf — -^fx e BsiSi te . (37) 
Here, we define the momentum transfer g = (go, <?) as 

g^Gfcl-lfc/l-AEfafci-fc/) (38) 
with AUq = U (f) — Uq(i), and rewrite the energy delta-function as 

5(N + ^(ft) + U (i) - \k f \ - E}( Pl + q)- U (f)) 

= S(E*( Pl ) + q - E}( Pl + q)) = 7-*-r5(t - t p ) , (39) 

Pi 9 



where i = q ■ pj{\q\ |p;|), and 



_ 2q E*{p i ) + q 2 + M?-Mf 

M\Pi\ ' ( j 

Furthermore we write 

aWte) + g - B;fe + «)) = - *p) + - f p) • ( 41 ) 

IPillyl PiQ 

Note that the terms proportional to s a (a = l,i,j) vanish in Eq. ( 128]) . and the 1^2,3,4 do 
not contribute to the final results to first order in fi b B. In view of this fact, we can further 
separate the magnetic part of the cross-section of Eq. (132)) into two parts as 

Aa = Aa M + Aa eU (42) 

11 



where the first and second terms are the contributions from the target particle and the 
outgoing electron, which appear only in the absorption (v e — > e~) process. 
The first term of Eq. ( I42p Ao^,/ is defined by 

dtdw k = i^rw { " n,(|fc/l)) a + b) (43) 



with 



T * = hi ^At-tv){<(E: + U^W-n f {E} + U {f))]^W 



where 



q\J \Pi\E* 

+ n' f (E} + U (f))n t (E* + U (t))( f x t W i - 2fi f W f ) 

Tb = ~ [ ^(Et + q )6(t-t p )n i (Et + U (i)) 
Q J PiEi 

x[l-n f (E} + U (f))] L^i-^^l) , (44) 
E*( Pt ) 

= c 2 v { [k f ■ (M* fPi - M?p f )] {h ■ h) - [h ■ (M* fPi - M*p f )] (k f ■ h)} 
+c\ { [-k f ■ (M} Pi + M*p f )] (hi ■ h) + [hi ■ (M* fPi + M*p f )] (k f ■ b t )} 
-2c v c A M* {(k f ■ p f )(ki ■ h) + (h-p f )(kf ■ bi)} , (45) 

JPfT + M f 

w f = v , , w f 

= c 2 v { [kf ■ (M* fPi - M*pf)} (h ■ b f ) - [h ■ (M* fPi - M*pf)} (k f ■ b f )} 
+c\ { [kf ■ (M) Vi + M* Pf )] (h ■ b f ) - [h ■ (Mj Pl + M*p f )] (k f ■ b f )} 
-2c v c A M} {(h ■ Pi )(kf ■ b f ) + (k f ■ Pi )(h ■ b f )} (46) 



with 



K = ^ P E>M" aaM ■ (47) 

In these equations the four momenta pt and pf are defined by pi = (E^p^p^ and pf = 
(E}(pf),p f ). 

When the target particle is an electron, the above expression is slightly altered. When 
both the initial and final particles are electrons, the above equations are written as 

Wi/m e = Cy {[kf ■ (p, t - pf)] (ki ■ bi) - [k { ■ (pi - pf)} (k f ■ bi)} 
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+c 2 A {[-k f ■ fa + p f )} (k ■ h) + [h ■ (pi + p f )) (kf ■ h)} 
-2c v c A {{k f ■ p f )(ki ■ hi) + (kf ■ bi)(ki ■ p f )} , (48) 
Wf/m e = c 2 v {[k f ■ (pi - p f )} (ki ■ bf) - [ki ■ (pi - Pf)} (k f ■ b f )} 
+4 {[kf ■ (P l + Pf)} (h ■ b f ) - [h ■ ( Pi + p f )} (k f ■ bf)} 
-2c v c A {(ki-pi)(kfb f ) + (kfPi)(ki-bf)} , (49) 

and 

Kf = T^-MPij)- (50) 
\Pij\ 

In the actual calculation we take the limit of m e — > 0, keeping n e Wi.f and bij finite. 
The second term in Eq. ( 1421 . Aa e i, is defined by 
-I -l 



4nG 2 F B 



(2 7 ry\q\\k i \\kf\ 

d 3 p 



d 3 . 
dkJ A(Xel 



« E n Klfc/l) / T^^At-tvm + ^WMEl + U^) [l-n f (E} + U (f))] 

+ ^li-njflfc,!)] / T^^-tpK^ + t/oWK^ + t/oa))^ 



[i- n/ (£; + w))]^ (si) 



with 



\kf\ 

-C V [(ki ■ Pf)( Pi ■ b e ) + (ki ■ Pi )(pf ■ b e ) - MfMi(h ■ be)) 
-c\ [(ki ■ Pf)(Pi ■ K) + (ki ■ Pi)(Pf ■ b e ) + MfMi(k ■ b e )} 

+2c v c A [(ki ■ p f )(pi ■ b e ) - (ki ■ Pi)(pf ■ b e )) . (52) 
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IV. RESULTS AND DISCUSSION OF NEUTRINO CROSS-SECTIONS 



In this section we present the cross-sections for neutrino scattering [y e — > v e ) and ab- 
sorption [y e — > e - ) in matter with and without a magnetic field. For the purposes of this 
investigation we set that the lepton fraction is fixed to be Yi = 0.4, and the neutrino incident 
energy is taken to be its chemical potential, \ki\ = e u , unless otherwise noted. We use the 
parameters of the weak-interaction, cy and ca, from Ref. 



A. Neutrino Cross-Sections without a Magnetic Field 
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FIG. 2: (Color online) Density dependence of the scattering (a and c) and absorption differential 
cross-sections (b and d) of neutrinos in neutron-star matter at T = 20 MeV without Lambdas (a 
and b) and with Lambdas (c and d). The initial momentum and angle of the incident neutrinos 
are taken to be = e v and 9i = 0°. Dotted, dashed and solid lines represent the results for ps = 
po, 3po and 5po, respectively. 

In Fig. |2] we show density dependence of the differential cross-sections per baryon for 
the scattering (a and c) and absorption (b and d) of neutrinos in matter without Lambdas 
(a and b) and with Lambdas (c and d) at T = 20 MeV. The subscripts 'S' and 'A' refer 
to the scattering or absorption cross-sections, respectively. Solid and dashed lines show the 
results in matter including Lambdas or no Lambdas, respectively. We see that the scattering 
cross-sections are forward peaked, while the absorption cross-sections decrease at forward 
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angles when pB < 3po- 
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FIG. 3: (Color online) Energy dependence of the scattering (a and c) and absorption differential 
cross-sections (b and d) when 0, = 0° for neutrinos in neutron-star matter at ps = 3/9o and T = 20 
MeV without Lambdas (a and b) and with Lambdas (c and d). Dotted, dashed and dot-dashed 
lines show the results for incident neutrino energies |fcj| = 100, 150, and 250 MeV, respectively. 
The solid lines show results for when the incident neutrino energy is set equal to the neutrino 
chemical potential, i.e. |fcj| = e v . 



In Fig. |3] we show the energy dependence of the differential cross-sections per baryon at 
Pb = 3po and T = 20 MeV for various incident neutrino energies. The solid lines show the 
results when the incident neutrino energies are set equal to the neutrino chemical potentials, 
i.e. \ki\ = e u = 196 and 217 MeV with and without Lambdas in matter. Dotted, dashed 
and dot-dashed lines represent the results at |fcj| =100, 150 and 250 MeV, respectively. 

When the incident energy is \ki\ = 100 MeV, the cross-sections show a minimum at 
forward angles. As incident energy increases, however, the cross-sections gradually become 
larger and finally become peaked at forward angles. This behavior arises from the the differ- 
ence in Fermi distributions between the spin-up and spin-down particles, as was discussed in 
Ref. [32j. This Pauli blocking affects the results at all angles, and, in particular, manifests 
itself at forward angles. However, this effect becomes smaller at higher incident energies 
as shown in Fig. El We have confirmed that the cross-sections always show forward peaks 
when we turn off the Pauli blocking term for the final lepton, (1 — n{). We can therefore 
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conclude that the effect of the Pauli blocking is clearly exhibited at low incident energy as 
a suppression of the differential cross-sections at forward angles. 



B. Differential Neutrino Cross-Sections in a Magnetic Field 

In this subsection we discuss effects of a magnetic field on the neutrino reactions in 
neutron-star matter. For illustration we first calculate the differential cross-sections per 
baryon, das t A/dQ/A with an initial neutrino angle of Q{ = 0° at a matter density of ps = 3po 
and a magnetic field of B = 2 x 10 17 G. This gives p^B = 0.63 MeV, where p^ is the nuclear 
magneton. The initial momentum is taken to be equal to the chemical potential in each 
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FIG. 4: (Color online) Effects of magnetic fields on the differential cross sections per baryon 
da/dil/A, from Eq. (|32p in units of 10 -16 fm 2 . This figure is for neutrino scattering (v e — > v e ) in 
neutron-star matter at a density of ps = 3/9o at T = 20 MeV (a) and 40 MeV (c). Lower panels 
(b and d) are the same as the upper panels but for neutrino absorption {y e — > e~ ). The initial 
momentum and angle of the incident neutrinos are taken to be \k{\ = e v and 6i = 0°. Solid and 
short-dashed lines represent the calculated results including Lambdas in the neutron star matter 
associated with and without a magnetic field of strength B = 2 x 10 17 G, respectively. Dot-dashed 
and dotted lines represent the results without Lambdas. 



In Fig. H]we show the neutrino scattering [y e — > u e ) cross-sections in the upper panels (a 
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and c) and the absorption {y e —> e~) cross-sections in lower panels (b and d), based upon 
Eq. fl32|) . The results in the left and right panels are for temperatures T = 20 and 40 MeV, 
respectively. Solid and dashed lines show the results for matter including Lambdas and no 
Lambdas, respectively. For reference, we also plot the results without a magnetic field both 
including (dot-dashed lines) and excluding (dotted lines) Lambdas. 
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FIG. 5: (Color online) Magnetic parts of the differential cross-sections per baryon dAa/dCl/A, 
which correspond to the 2nd term in Eq. (|32p . in units of 10 -16 fm 2 for neutrino scattering {y e — > u e ) 
on neutron-star matter at ps = 3po for T = 20 MeV in a system without Lambdas (a) and with 
Lambdas (c). Lower panels (b and d) are the same as the upper panels but for neutrino absorption 
{y e — > e~) in a system without Lambdas (b) and with Lambdas (d). The initial momentum and 
angle of the incident neutrinos are taken to be = e u . Solid, dashed and dotted lines represent 
the calculated results when 6{ = 0°, 90° and 180°, respectively. 

This figure indicates that the magnetic field does not affect much the scattering cross- 
sections when B w 2 x 10 17 G. The contribution from each individual particle such as protons 
or neutrons is not so small. These contributions, however, tend to cancel each other. Nev- 
ertheless, the magnetic field suppresses the absorption cross-section in the forward direction 
and enhances it in the backward direction. In particular, near 9f ~ 0°, the suppression from 
the magnetic field is as much as 20 —30 %. This contribution is almost as large as that from 
the A particle which has strangeness. 
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In Fig. [5] we show the magnetic parts of the differential cross-sections, Act of Eq. (|32|) . at 
a matter density of ps = 3po and a temperature of T = 20 MeV. The upper panels are for 
neutrino scattering [y e — > v e ) and the lower panels (b and d) are for absorption [y e — > e - ). 
The results in the right and left panels are for matter including and excluding Lambdas, 
respectively. Solid, dashed and dotted lines represent the calculated results for incident 
angles, of Q\ = 0°, 90° and 180°, respectively. In these calculations we keep the difference of 
the azimuthal angle between the initial and final leptons equal to be zero, i.e. 0/ — <pi = 0. 

These calculations show that the magnetic field enhances the scattering cross-sections in 
the direction along the magnetic field (arctic direction). For absorption an enhancement 
is seen in the opposite direction (antarctic direction). This role of the magnetic field in 
making an asymmetry of the scattering and absorption cross-sections of neutrinos inside 
neutron-star matter would lead to the coherent effect of enhancing the neutrino drift in the 
arctic direction while suppressing it in the antarctic direction, as will be discussed below. 

C. Angular-integrated Neutrino Cross-Sections in a Magnetic Field and The 
Asymmetries 

In order to discuss the effects of neutrino transfer inside the PNS, we here calculate the 
scattering cross-sections integrated over the momenta of the initial neutrinos as 



The absorption cross-sections are however integrated over the momenta of the final electrons 

as 



Note that the non-magnetic parts of the integrated cross-sections, A , also follow the same 
equations. 

Figures [6] and [7| show Aas/cr^ with |fej| = e u and Aff^/aJJ with \kf\ = e u as functions 
of 8f and 6^, respectively, for various matter densities in the interval, p < Pb < 5p - We 
plot the calculated results for matter without Lambdas (upper panels) and with Lambdas 
(lower panels) at T = 20 MeV (left panels) and T = 40 MeV (right panels). Similar to 
the differential cross-sections, the magnetic field enhances the neutrino integrated scattering 
cross-sections and suppresses the integrated absorption cross-sections in the arctic direction 




(53) 




(54) 
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FIG. 6: (Color online) Ratios of the magnetic part of the scattering cross-sections Aug, to the cross- 
sections without a magnetic-field on matter without Lambdas (a and c) and with Lambdas (b 
and d) at T = 20 MeV (a and b) and at T = 40 MeV (c and d). Dotted, dashed and solid lines 
represent the results for densities of ps = Po, 3po and 5po> respectively. 
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FIG. 7: (Color online) Ratios of the magnetic part of the absorption cross-sections Acta, divided 
by the cross-sections without a magnetic-field a\ on the matter without Lambdas (a and c) and on 
matter with Lambdas (b and d) at T = 20 MeV (a and b) and at T = 40 MeV (c and d). Dotted, 
dashed and solid lines represent the results for ps = po, 3po and 5po, respectively. 
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parallel to the magnetic field B. The magnetic field has an opposite effect in the anti- 
parallel antarctic direction. This increases the neutrino emission in the arctic direction and 
decreases it in the antarctic direction. 
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FIG. 8: (Color online) Contributions from each particle process to the magnetic part of the scatter- 
ing cross-sections at T = 20 MeV divided by the integrated cross-sections without a magnetic-field, 
<jg. Left upper (a), left lower (b), right upper (c) and right lower panels (d) exhibit the results 
without and with Lambdas at a density of ps = 3po and pb = 5po, respectively. Dashed, solid, 
dot-dashed and long dashed lines represent contributions from electrons, neutrons, protons and 
Lambdas, respectively. Dotted lines represents sums of the contributions. In the right bottom 
panel (d) the solid and long dashed lines cannot be distinguished. 

In Fig. [8], we show each contribution to the scattering cross-sections from the constituent 
particles of matter without Lambdas (upper panels) and with Lambdas (lower panels) at the 
densities of ps = 3po (left panels) and ps = 5po (right panels). The contribution from the 
protons is opposite to those from electrons, neutrons and Lambdas because of the different 
signs of the magnetic moments. These contributions tend to cancel to each other, and the 
magnetic parts of the scattering cross-sections become slightly smaller. However, when one 
allows Lambdas to appear in the system, the proton fraction decreases and in this case the 
cancellation is not as large as the case without Lambdas (see Fig. [[]). 
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FIG. 9: (Color online) Contributions from each particle process in the magnetic part of the ab- 
sorption cross-sections in matter with Lambdas at T = 20 MeV at ps = 3po in upper panels and 
at pb = 5po in lower panels. Left and right panels exhibit results divided by the non-magnetic 
parts of the total cross-sections and their non-divided respective contributions. Solid and dashed 
lines represent the contributions from the n — > p and A — > p processes, respectively. 

The solid and dashed lines in Fig. [9] show the contributions from the n — > p and A — > p 
neutrino absorption processes, respectively. Upper and lower panels exhibit the results at 
Pb = 3p and Pb = 5p , respectively. Results in the left and right panels are divided by 
the non-magnetic parts of the integrated cross-sections and their respective non-magnetic 
contributions. 

Contributions from the A — > p process seem to be much smaller than those from the 
n — > p in the left panels, but in the right panels the former contributions are as large as 
the latter. This apparent difference is because the Cabibbo angle is very small, sin 2 9c ~ 
5.0 x lCT 2 . Namely, because the non-magnetic part of the A — y p process is associated with 
the strangeness change (AS = 1), its transition probability is ~ sin 2 9 C times smaller than 
that of the n — >■ p process with no strangeness charge (AS* = 0). As a result, contributions 
from the A — > p process to the total non-magnetic part becomes very small. However, 
when one divides the small contributions from the magnetic part of this process by the 
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small quantities of those from respective non-magnetic parts, the ratio shows an interesting 
difference as illustrated in the right panels. With Lambdas present in the matter, the proton 
fraction becomes smaller as the density changes, and the contribution from the magnetic 
parts of the A — > p process becomes remarkably large. 

D. Neutrino Mean- Free- Paths 
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FIG. 10: (Color online) Upper panels (a) and (c) show the neutrino MFP for scattering (a) and 
absorption (c) in neutron-star matter without a magnetic field. The lower panels (b) and (d) show 
the magnetic contribution to the MFP for scattering (b) and absorption (d). Since the magnetic 
part of the MFP for scattering is negative, we multiply by (—1). Solid and dashed lines represent 
the results in matter at T = 20 MeV both with and without Lambdas, respectively. Dot-dashed 
and dotted lines represent the results in matter at T = 40 MeV with and without Lambdas, 
respectively. 

In order to consider the application of the above results to astronomical phenomena, we 
discuss the neutrino mean- free-paths (MFPs). In Fig. [TO] we show the density dependence of 
the neutrino MFPs Xs,a = V/ as,A for the scattering (a) and the absorption (c) processes for 
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B = at T = 20 and 40 MeV in the upper panels, where V is the volume of the system. For 
this illustration, the incident neutrino energy is fixed to be equal to its chemical potential. 

The scattering and absorption MFPs rapidly decrease as the density increases up to 
Pb ~ (2 — 3)po- When the system does not include Lambdas, both MFPs (dashed and 
dotted lines) decrease monotonically even beyond p# ~ (2 — 3)p . When the system includes 
Lambdas, the scattering mean-free-paths also decrease, but the absorption mean-free-paths 
increase in ps > 3po, because the cross-sections for p e + A — > p + e~ are smaller than those 
of v e + n — > p + e~ . 

In addition, we show the magnetic 
contributions to the MFPs, A\s,a = 
[V/a(0°) - V/<7(180°)]/2, in the lower pan- 
els (b and d). Since As < 0, we multiply 
by (—1). We see, again, that the contribu- 
tion of the magnetic field is ~ 1 — 2% of 
the non-magnetic parts. 

The slopes of the magnetic parts of neu- 
trino cross-sections A<js,a are almost con- 
stant as function of cos O^j (see Figs. [6] 
and [7]). Hence if we define the slopes as 

S s ,a = {A- (7 s,a/o's,a)/ cos ®i,fi Acr 5, then 
the integrated cross-sections Aas,A can be 
approximately written as 
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FIG. 11: (Color online) Density dependence of Ss 
when |fej| = e v (a) and Sa when \kf\ = e u (b). 
The discrepancy in the use of this formula Dashed, dotted, solid and dot-dashed lines rep- 
is less than 1 %. resent the results in matter without Lambdas at 
Since S s > and S A < 0, the neutrinos T = 20 MeV, without Lambdas at T = 40 MeV, 
scatter strongly and absorb weakly in the w ith Lambdas at T = 20 MeV and with Lambdas 
arctic direction due to the magnetic field. a t T = 40 MeV. 
In Fig. CD], we show the density dependence 

of Ss (a) and Sa (b). It is evident that the effects of the magnetic field become smaller as the 
temperature and the density increase. This density dependence arises from the fact that Act 
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is approximately proportional to the fractional area of the distorted Fermi surface caused 
by the magnetic field. Hence, the relative strength A<js,a/ &s a diminishes with increasing 
density. 

However, the density dependence of S$ in the system including Lambdas exhibits a local 
minimum around p# ~ 3 and increases again in the density region, 3p ^ Pb ^ (5 — 6)p . As 
commented before, the Lambda fraction rapidly increases for p# > po, and its contribution 
enhances Ss (see Fig. [I] and Fig. [8]). 

E. Neutrino Transport and Pulsar Kick Velocities 

In this subsection we discuss implications of these findings for neutrino transport in 
strongly magnetized PNS. We estimate the momentum transfer from the asymmetric neu- 
trino emission. 

1. Boltzmann Equation 

To estimate the neutrino momentum transport we begin with the phase-space distribution 
function f u (r,k) and calculate the asymmetric neutrino emission. This f v satisfies the 
following Boltzmann equation: 

(l + L ^)^ (r ' fc)=/co " (56) 

with 

- fMhiPi) [1 - fi{h)\ [1 - / 2 (P 2 )]} , (57) 

where Wif is the reaction probability. The index / denotes lepton, electron or neutrino, and 
the indices 1 and 2 label the target particles, e.g. baryons and electrons. 

Here, we introduce several assumptions to obtain a solution to the Boltzmann equation 
as follows. First, we assume that the system is almost in equilibrium, and that f u {f, k) can 
be separated into two parts 

Mr, ft) = Mr, ft) + Af(r, k) = - ' - + Af(r, k) , (58) 
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where the first and the second terms are the local equilibrium part and the deviation from 
the equilibrium, respectively, with e u (r) the neutrino chemical potential at the position 
r. The phase-space distribution functions for the other particles are assumed to be local 
thermodynamic equilibrium distributions. In addition, we also omit the contribution from 
e~ + B — > B' + v e . The collision term can thus be written as 

r°* - / (w s {A f ( kl ) [f 2 (i - / (fc))(i - /i) + / (fe)A(i - h)\ 

- Af(k) [/o(feO/ 2 (l - fx) - A(l - /o(*i))(l - / 2 )]} 
-W A A/(fc) [/!(!-/.(*,))(! "/a)]) , (59) 

where W5 and Wa are the scattering and absorption probabilities. 

We make further assumptions that only the absorption process makes a dominant con- 
tribution to the neutrino momentum transport. This is because the magnetic contribution 
to the scattering cross-section is small. 

Under these assumptions, and when the system is in steady state, df u /dt = 0, the 
Boltzmann equation for the neutrino transport can be expressed as 

where the absorption cross-section, a a, is a function of k and Pb{ t )- 

In the present approximation the neutrinos are taken to propagate along a straight line, 
and the above Boltzmann equation can be analytically solved as follows. First, we define 
a plane A Q that is perpendicular to the neutrino momentum k. This plane is constructed 
to intersect the center of the neutron star, which we take to be the origin of the coordinate 
system r = (0, 0, 0). Then, we introduce xl and Rt such that r = + Rt, where xl is 
the component of r parallel to k and Rt -L k. In terms of xl and Rt, Eq. (1601) can then 
be written as 

de v df dAf °A Aft D . , (MS 

dx- L dT„ + ^; = -v Af ^ RT > k) ' (61) 

where Rt = \Rt\ and de u /dxL = (k ■ f)de u /dr. The solution is given by 



r x L [ de df 1 I f XL a a 

A f (x L ,R T ,k) = j o d y['didr u r p [-J y dz v 



(62) 
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As neutrinos are created inside a PNS and 
propagate through the matter, their intensity 
will be attenuated by absorption. The expo- 
nential in Eq. f[6"2"j) accounts for this feature. 
For sufficiently large <Ja/V, we expect that 
very few neutrinos produced deep inside the 
PNS will reach the surface. 

subsubsectionMean free path in NS matter 
To give a more concrete picture we next ana- 
lyze the MFP of neutrinos. Fig. [TJ] shows the 
neutrino MFP for absorption = (cja/V)~ 1 
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path for absorption with a neutrino energy E u = 
30 MeV (solid line), 50MeV (dashed line) and 
lOOMeV (dotted line) in neutron-star matter at 
T = 20 MeV without a magnetic field. The thick 
and thin lines represent the results in matter 
with and without Lambdas, respectively. 



(dotted line) in neutron-star matter at a tem- 
perature of T = 20 MeV without a mag- 
netic field. The thick and thin lines rep- 
resent the results in matter with and with- 
out Lambdas, respectively. The MFPs are 
smaller than a few km so that most of the 
neutrinos produced in the central region are 

absorbed. However, the neutrinos produced near the surface contribute to the net emission 
of neutrinos; this fact is qualitatively the same as the result obtained in Ref. _5(3]. In addi- 
tion, we see that the neutrino MFP is longer when its energy is large because of the Pauli 
blocking effect of the final electron. As a result lower energy neutrinos are absorbed more 
efficiently. 

In order to solve Eq. fl62l) . we need to know a^/V as a function of the density p, the 
magnetic field B, the initial neutrino energy E v , and the angle between the magnetic field 
and the initial neutrino momentum, 9 V . For this calculation we have made a data base of 
(J A ctS £1 function of the baryon density ps and the incident neutrino energy E v . 

However, it is not easy to make a data base of the magnetic part of A<tm because it is 
a function of pb, E v and 9 V as well as B. This leads to a computationally intensive five 
dimensional integration. Therefore, we introduce a fitting function for the magnetic part 



26 



deduced as follows. 

As defined by Eq. fl55|) . the angular dependence can be approximately written as 

a A = a A (l + S A cos9„). (63) 
We then assume that S A obeys the following function: 

-S A = A M + B M e- CMEl (64) 

with 

A M = A + A, O^Y , (65) 

The values of parameters are given in Table I. All quantities except ps and E v are constant 
and adjusted to reproduce the theoretical results shown in Fig. [13] as described in the figure 
caption. 



2. Proto Neutron-Star Model 



Baryon density profiles of our PNS model are shown in Fig. [TH We assume an isothermal 
PNS model at temperatures T = 20 and 30 MeV with and without Lambda particles. 
For this illustration we fix the total gravitational mass of the PNS to be 1.68 M & . The 
appearance of Lambda particles in ps > 2p softens the EOS. This increases the baryon 
density and the neutrino chemical potential. The density profiles with Lambdas are sensitive 
to the temperature. 



3. Momentum transfer 

We use these density distributions of the PNS to calculation the neutrino momentum 
transport. We define the effective spherical surface Sn where ps = Po, and estimate the 
kick velocity from the angular dependence of the emitted neutrino momentum at this surface. 
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3.92 x 10" 1 


1.57 x 10~ 2 


8.75 x 10~ 2 
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0.355 


0.3920 


0.354 


0.0972 
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-2.96 x 1(T 3 


2.62 x 10~ 3 


-5.35 x 10- 3 


-2.01 x 10~ 3 


Bi 


-1.23 x 1(T 3 


-4.02 x 10~ 3 


5.15 x 10" 4 


-9.81 x 10~ 4 


B 2 


-1.72 x 10" 4 


2.03 x 10~ 4 


-1.89 x 10" 4 


2.36 x 10" 5 


B3 


-8.43 x 1(T 7 


-9.96 x 10~ 7 


-9.26 x 10~ 7 


1.16 x 10~ 7 


C (MeV- 2 ) 


1.16 x 1(T 5 


1.05 x 10~ 5 


2.51 x 10~ 5 


4.31 x 10" 7 


Ci(MeV- 2 ) 


-3.90 x 10~ 8 


4.37 xl0~ 7 


-7.17 x 10~ 6 


4.99 x 10~ 7 


C 2 (MeV- 2 ) 


1.63 x 10~ 7 


9.67 x 10~ 8 


7.46 xl0~ 7 


4.30 xl0~ 9 


C 3 (MeV- 2 ) 


-5.59 x 10~ 6 


-4.23 x 10~ 9 


-2.28 x 10~ 8 


-4.42 x 1Q- 10 



TABLE I: Parameters of Eqs. (1651) - (16 



The total momentum per unit time of the neutrinos emitted along the direction n is then 
calculated as 



P 



dr 



d 3 k 



Af(r,k)(k-n)S(k- (fc-n)n) 



>s N J (2tt,- 
The momentum P can be approximately written as 

P = P + AP ^ P + P! cosd 



(68) 



(69) 



in terms of the polar angle 9. The asymmetry of the neutrino momentum, AP/P , as a 
function of 6 is shown in Fig. [T51 

We use the distribution in (169]) to calculate the ratio between the average momentum in 
the direction of the magnetic field < P z > and the total emitted neutrino energy E r p 1 i.e. 
< P z > /E T = Pi/3P . The calculated results are < P z > / E T = 0.0203 and 0.0196 in 
systems with and without the Lambdas at T = 20 MeV, and < P z > /Ex = 0.0078 and 
0.0092 in a system with and without Lambdas at T = 30 MeV, respectively. 



We assume that the total energy emitted in neutrinos is Ex ~ 3 x 10 53 erg 3jJ. For the 
Mns — 1-68 M Q isothermal model with T = 20 MeV, the calculated kick velocities tare 
Vkick =< Pz > /Mns — 610 km s~ x or 580 km s~ 4 in neutron-stars including Lambdas or no 
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FIG. 13: (Color online) The magnetic slope parameter of the neutrino absorption versus incident 
energy at T = 20 MeV. Open and full circles show the results in the present calculation at ps = Po 
in the system without and with Lambdas, respectively. Open and full squares indicate those at 
Pb = 3po hi a system without and with Lambdas. Dotted, dot-dashed, dashed and solid lines show 
the results of the fitting function in the system at pb = Po without and with Lambdas, and those 
at pb = epo without and with Lambdas, respectively. 



Lambdas, respectively. At T = 30 MeV, their velocities are Vkick =< P z > /M NS = 230 km 
s _1 and 270 km s _1 . 

In actual observations, the average value of the kick velocity is about v^ick — 400 km 
s _1 , and the highest reported value is ~ 1500 km s _1 . Our calculated values are thus close 
to the observed average pulsar kick velocity. We note that Lai and Qian [31] obtained a 
similar result (vkick =280 km However their result was calculated in a non-relativistic 
framework and without Lambda particles. 

The existence of the Lambdas has a small but non-negligible effect on the pulsar velocities. 
In Fig. [TT]we can also see that the asymmetry of the absorption cross-section becomes smaller 
when the temperature increases. This is because the smearing of the distribution functions 
causes the n' b term in Eq. f )44|) to be small. Hence, the kick velocities at T = 30 MeV are 
smaller than those at T = 20 MeV. 

In the central region, high energy neutrinos E v > 100 MeV are copiously produced, but 
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FIG. 14: (Color online) Density distribution versus the PNS radius. Solid and dashed lines repre- 
sent the results in matter with and without Lambdas at T = 20 MeV, and dot-dashed and dotted 
lines represent the results in matter with and without Lambdas at T = 30 MeV, respectively. 
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FIG. 15: (Color online) The variation of emitted neutrino momentum versus the polar direction. 
Solid, dot-dashed, dashed and dotted lines represent the results in a system with Lambdas at 
T = 20 MeV, and T = 30 MeV, and those without Lambdas at T = 20 MeV and T = 30 MeV. 



their MFP is only about several 10 3 cm. They are, therefore, almost completely absorbed in 
the transport process. Most of neutrinos with energy < 50 MeV contribute to the momentum 
kick because the MFP for these neutrinos is larger (Fig. [T2l) . If one can assume that the 
thermalization process is faster than the time scale at which the neutrino absorption process 
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directly affects the collective motion of the PNS, then only information on the cross-section 
in the low density region affects the final asymmetry of the neutrino emission. 

Neutrinos are continuously further absorbed in the lower density regions before they are 
emitted outside the neutron-star, and the asymmetry is should be retained. Indeed, when we 
extend the calculation to much lower density ps = 0.5p , we find that the the asymmetries 
are almost the same as the above results, but that the energy of the emitted neutrinos is 
small. 

We caution, however, that in such low density regions, both the magnetic field and 
temperature may be lower than those assumed in the present isothermal model. If instead 
of an isothermal neutron-star model, one were to use a an isentropic model with uniform 
entropy, then the kick velocity may be smaller. In the surface region of magnetars, the 
magnetic field is still as high as 10 15 G, however that is only about 1/100 of the value 
adopted in the present calculation. A lower magnetic field may reduce ifedb but the lower 
density and temperature may enhance it. In such a subtle situation the scattering process 
which we ignored in the present calculation should also be included as it enhances the 
neutrino asymmetry. This could tend to increase the kick velocity. 
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V. SUMMARY AND DISCUSSION 



We have studied the neutrino scattering and absorption processes in strongly magnetized 
proto-neutron stars (PNSs) at finite temperature and density. We used a fully relativistic 
mean field (RMF) theory for the hadronic sector of the EOS including hyperons. We solved 
the Dirac equations for all constituent particles, p, n, A, e, and u, including a first order 
perturbation treatment of a poloidal magnetic field with B ~ 10 17 G. We then applied 
the solutions to obtain a quantitative estimate of the asymmetry that emerges from the 
neutrino-baryon collision processes. We took into account the Fermi motion of baryons 
and electrons, the momentum dependence of their spin vectors, their recoil effects, and the 
associated energy difference of the mean fields between the initial and final particles exactly. 
We thus included the most important effects of the distortion of the Fermi spheres made 
by the magnetic field in this fully microscopic framework that implies asymmetric neutrino 
scattering and absorption cross-sections. 

We found that the differential neutrino absorption cross-sections are suppressed in the 
arctic direction parallel to the poloidal magnetic field B in both cases with and without 
As, while the differential scattering cross-sections are slightly enhanced. On the other hand, 
as expected from the sign of the couplings between the magnetic moments of baryons and 
the external field, the neutrino absorption and scattering cross-sections are respectively 
enhanced and suppressed in the antarctic direction. This is completely opposite to those in 
the arctic direction. The differential cross-sections were integrated over the momenta of the 
final electrons for absorption and over the momenta of initial neutrinos for the scattering, 
respectively. Quantitatively, when B = 2 x 10 17 G, the reduction for the absorption process 
is about 2%, and the enhancement for the scattering process is about 1% in the forward 
direction along the direction of B. 

Several interesting facts are evident in the angular distributions of both cross-sections, 
which depend on the magnetic field, the baryon density, and the temperature of the PNS 
matter. Among them, we find, an appreciable forward suppression and backward enhance- 
ment in the differential absorption cross-sections due to the difference in Fermi distributions 
between the spin-up and spin-down particles. This effect is larger at lower neutrino incident 
energy. The asymmetry becomes smaller as the density increases because the asymmetry 
arises from the magnetic part of the cross-sections which is proportional to the distortion 
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of the Fermi surfaces caused by the magnetic field. This tends to diminish with increasing 
matter density. 

Using these cross-sections, we calculated the neutrino mean-free-paths (MFPs) as a func- 
tion of the baryon density and temperature within a PNS. We then applied the above results 
to a calculation of pulsar-kicks in core-collapse supernovae. We solved the Boltzmann equa- 
tion using a one-dimensional attenuation method, assuming that the neutrinos propagate 
along an approximately straight line and that the system is in steady state. We only in- 
cluded the MFPs for neutrino absorption which dominates over scattering in producing the 
asymmetric momentum transfer to the PNS. 

We estimated pulsar kick velocities from the calculated total momentum per unit time 
that is transferred from the emitted neutrinos to the PNS along the direction parallel to 
the poloidal magnetic field B. For an isothermal neutron-star with M^s = 1.68M Q and 
a total energy in emitted neutrinos of E T 3 x 10 53 erg, the estimated kick velocities are 
v kick = 610 km s' 1 and 580 km s' 1 (at T = 20 MeV) or v kick = 230 km s" 1 and 270 km 
s _1 (at T = 30 MeV) including Lambdas or no Lambdas, respectively, These values are in 
reasonable agreement with the observed average pulsar-kick velocity of Vkick = 400 km s _1 . 

VI. FUTURE WORK 

In the present calculations we have adopted several assumptions which we summarize 
here both as a caveat for the reader and as a summary of issues to be addressed in future 
work. One such assumption was ignoring the neutrino scattering process. This scatter- 
ing might enlarge the kick velocity. The one-dimensional attenuation method to solve the 
Boltzmann equation is also a coarse approximation. We have assumed that the asymmetry 
in neutrino emission is dominated by the emission from low-density regions with ps < 3po 
where neutrino opacity changes drastically. And, we have assumed that the internal high- 
density region only contributes to the neutrino diffusive flux which smears out the expected 
neutrino asymmetry. However, as was discussed in the last section, if the thermalization 
process is considered dynamically the asymmetric neutrino scattering and absorption in the 
high-density region might also contribute to an aligned drift flux along the direction of B. 
This could generate a gradual acceleration of the pulsar-kick. Numerical simulations of the 
neutrino transport inside the PNS coupled to our microscopic calculations of the asymmetric 
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neutrino scattering and absorption cross-ksections are highly desirable in order to address 
to this critical question. It has been pointed out by Arras and Lai 3(3] that the neutrino 
distribution tends to be asymmetric only near the surface of PNS. This is consistent with the 
picture adopted in the present attenuation approximation for the neutrino transport. The 
issue becomes more subtle, however, if the thermalization process is considered. It would be 
interesting to clarify by numerical calculations the extent to which the asymmetric neutrino 
scattering and absorption contribute to the drift velocity as well as the diffusive velocity 
of outgoing neutrinos considered here. Another important question is to address the link 
among asymmetric neutrino-baryon collisions, neutrino drift, and the collective response of 
the PNS as pulsar kick. Further investigations must be done by numerically solving the 
Boltzmann equation for the neutrino transport inside a PNS without approximations, al- 
though we believe that our adopted scheme of attenuation is more or less consistent with 
the microscopic picture. Numerical calculations including several dynamical effects are now 
underway. 

We also should take account of neutrino reactions in the much lower density region, 
Pb 3> Po, although we did not include that in the present study because of the numerical 
difficulty in calculating thermodynamic quantities of the EOS in the RMF theory. In such 
low density regions, the magnetic field is weaker, but the width of the Landau level, \/2eB, 
could be of the same order as the electron Fermi momentum, and it may affect the neutrino 
reactions. 

The strength of the magnetic field inside the PNS can easily reach 3 — 4 x 10 18 G in the 
high-density region according to the scalar virial theorem. This could make considerable 
effects, and a non-perturbative treatment of the magnetic field must be applied for this high 



field strength 



39]. We may again need to take account of the Landau levels. 



We also did not take account of the resonant spin-flavor conversion 44( in the magnetized 



451 ] or the self-interaction 



PNSs, and the neutrino-flavor conversion due to the MSW effect 
effect 46( in the present calculations. All of these could alter the asymmetric neutrino 



emission. A quark-hadron phase transition |47| or a hyper-nuclear matter phase [48] un- 
der a strong magnetic field is also considered to be another source to affect the neutrino 
asymmetry. 

If a poloidal magnetic field exists in the progenitor stars for SNe, a stable toroidal mag- 
netic field also is created in the core-collapse and explosion. In this case, the angular de- 
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pendence of the neutrino reactions may show more complicated interesting behavior. Thus, 
there are many open questions to be addressed in the future studies which are beyond the 
scope of the present article. 
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Appendix A: Dirac Spinor in a Magnetic Field 



In this appendix, we explain the detailed expressions of the Dirac spinor under the in- 
fluence of a magnetic field. The Dirac spinor u(p) can be obtained by solving the following 
Dirac equation 

K(p)u(p, s) = \i-M- U (b) - U T a z ] u(p, s) = 0, (Al) 
where Ut = fiB. Since the Green function S(p) is defined as 

K(p)S(p) = l , (A2) 

this Green function is written as 

S(p) = det K{S + S 1 U T + S 2 U^ + 5 3 C/|) , (A3) 

with 

det K = Pq — 2p 2 (p 2 + M 2 + E/|) + (p 2 + M 2 ) 2 + 2U^(p 2 z - p 2 T - M 2 ) + U$ , 

5 = (p 2 -E 2 )(i> + M) , 

51 = {j>l + El)a z + 2Mp a zlQ -2p z {j}-cT) lQ 

+ 2Mj9 2 7 5 7o + 2ip Pj/7°7 1 - 2ip Q p x ^rf , 

5 2 = -po7° + P 2 7 3 - PX7 1 - Pyl 2 + M , 

53 = -a z . (A4) 

Here the single particle energy of this Dirac spinor is obtained from det K = and 
becomes 



e(p, s) = ^pl+(^/^+AP + sU T ^ = sJe* + 2sU t ^ 2 t + M2 + U 2 , (A5) 

where s — ±1, and E p = a/p 2 + M 2 . Then, det K is rewritten as 

detK = (p 2 - e 2 (p, l))(p 2 - e 2 (p, -1)) . (A6) 
Furthermore, the Green function for this particle is written as 

Ujp, S)u(p, S) y(-p 

±1 Po - e(p, s) ± i<J ^ p + e(p, s) + i«5 



sip) = k-\p) = y: u(p, 1 Hp i\, + y, v }~ p ^ { : p :.] , (a?) 
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where u(p, s) and v(—p, s) are the Dirac spinors of the positive and negative energy states, 
respectively. 

By using the above quantities, we can obtain the Dirac spinor as 



u{p, s)u{p, s) = lim (po - e(p, s))S(p) 

P0->e(p,s) 



(A8) 



Now we expand S with respect to Ut and determine the Dirac spinor in first order 
perturbation theory. Here we define 

Po - e(p, s) 1 



D P 



lim 



(A9) 



P0 ^e( P ,s) detK 8e(p,s) {sU T y/Pr + M 2 

When \Ut\ *C 1, we can substitute p = e{p, s)&E p + sU T ^p 2 T + M 2 /E p into Eq. flAi]) 



and obtain 

D e S 



s_( 1 _ sU T ^/p 2 T + M 2 \ ( 



sU 7 



4£„ 



E 2 

p 



2^p 2 T + M 2 



... trs sU T Jp 2 + M 2 
x ^ (jj + M) + TVj 7 7o 



4£„ 



(p 7 + M) + 



-En 



7o 



,2 ^2 



st/ 7 



(AlO) 



Po=E p 



U T D e S x 



SE p ^p 2 T + M 2 



i sU T ^p 2 T + M 2 



El 



. . sJpl + M 2 
x <Sx + 2{E p a z + M<y cr z + ip y a x - ip x a y )— — U T 



E n 



p =E p 



1 



Si 



4E p [ + M 2 

+ U T 

1 



S 1 

' + —(E p a z + M^o<?z + ip*^ - iPy<7y) 



AE„ 



2E 2 E p 

Vz 



p =E p 



U T D e S 2 



s(i> + M) 75 ^ + ^(/3(r-p- M 75 ) U T 
v 



8E p ^/p 2 T + M 2 



Po=E p 

(-E pl ° + M + p z7 3 - p^ 1 - p yl 2 ) 



with 



1 



yjp\ + M 2 



( Pz ,0,0,E p 



(All) 



(A12) 
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Then, the Dirac spinor is written up to the first order in Ut as 

x / x + M)(l + i 5 4(p)s) p z U T . 
u(p,s)u(p,s) « ^ 4jgp + IgF (°- • P ~ M ^7o) 



' (-£ p7 ° + M + p 27 3 - - Pyl 2 ) . (A13) 
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